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Abstract 

The phase structure of the bosonized multi-flavor Schwinger model is investigated by 
means of the differential renormalization group (RG) method. In the limit of small fermion 
mass the linearized RG flow is sufficient to determine the low-energy behavior of the A^- 
flavor model, if it has been rotated by a suitable rotation in the internal space. For large 
fermion mass, the exact RG flow has been solved numerically. The low-energy behavior of 
the multi-flavor model is rather different depending on whether = 1 or > 1, where N 
is the number of flavors. For N > 1 the reflection symmetry always suffers breakdown in 
both the weak and strong coupling regimes, in contrary to the = 1 case, where it remains 
unbroken in the strong coupling phase. 

Key words: Renormalization group; Field theories in dimensions other than four 
PACS: ll.lO.Gh, ll.lO.Hi, ll.lO.Kk 
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Introduction. — Two-dimensional quantum electrodynamics (QED2) or the 
Schwinger model [1] exhibits many analogies with four-dimensional quantum 
chromodynamics (QCD4) including confinement, chiral condensate, topological 9- 
vacuum. The Lagrangian of QED2 with massive iV-flavor fermions which is called 
the iV- flavor (or multi- flavor) Schwinger model reads [Q,[3l,3,|5|] 



N 



1 



n=l ^ 



(1) 



where F^i, = d^A^, — duA^. Using bosonization technique [Q,[3l,Bj5ij6 Jvjsithe 
fermionic theory dB can be rnapped onto an equivalent Bose form 



(2) 
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with jj?' = 1 71, c = e"' / (27r) where 7 = 0.5774 is the Euler's constant, 9 is the vac- 
uum angle parameter, A^^ denotes normal-ordering w.r.t. m and y9„ n = 1, are 
one-component scalar fields. Both the fermionic and the bosonic form of the model 
has been analyzed by various methods from various aspects, e.g. mass perturbation 
theory [12], density matrix renormalization group (RG) method [fioll . lattice cal- 
culations [lid. Il4l ll5n. momentum RG method i20ll . etc. Physical properties (like 
e.g. chiral condensate [[s, 13, 15, 16, 17], boson mass spectrum [14,13]) have been 
investigated for arbitrary values of 9, fermion mass and temperature. 

The N = 1-flavor Schwinger model for 9 = ±tt has two phases [HSEIulQ. 
Illustrative and detailed analysis of the phase structure is presented in [TT]. The 
behavior of the theory is controlled by dimensionless ratio g/m. For g/m large, i.e. 
for strong coupling, the symmetry ^ —ip is unbroken, there is a unique vacuum 
at 99 = and there are no half-asymptotic particles. For g/m small, i.e. for weak 
coupling, the reflection symmetry suffers spontaneous breakdown, there are two 
vacua approximately located at (/? = ±y^/2 and half- asymptotic particles appear. 



The multi-flavor (A^ > 2) model has not been studied as extensively as the 1- 
flavor model. However, the relative ignorance toward the multi-flavor Schwinger 
model is perhaps not fully justified as it shows more resemblance to QCD4, because 
the model features a chiral symmetry breakdown. Based on the study of chiral 
condensate ['liflllllTilisI], the behavior of the Schwinger model has been found to 
be distinctively different for = 1 and for N > 2. Recently, the phase structure 
of the A^ = 1 flavor Schwinger model has been investigated by exact functional 



RG method [|21|] and it has been recovered in a rather straightforward way. Our aim 
in this work is to extend the RG analysis for the bosonized multi-flavor (A^ > 2) 
Schwinger model to consider its phase structure and clarify the difference between 
the 1-flavor and the multi-flavor models. 



The multi-flavor Schwinger model has relevance in solid state physics, too. It has 
been used to describe antiferromagnetic spin chain (see e.g. and the Bose form 
of the 1-flavor model has been proposed as an adequate model for the description of 
the vortex properties of two-dimensional (2D) isolated thin superconducting films 
[|23|] and the multi-flavor model has been used for description of vortex dynamics in 
magnetically coupled layered superconductors [|24|] . The number of flavors of the 
bosonized multi-flavor Schwinger model is equal to the number of layers of the su- 
perconducting layered system and the Fourier amplitude of the bosonic model ^ 
is related to the fugacity of the vortex gas. The RG techniques, like the real space 
RG method developed for spin systems usually rely on the dilute gas approximation 
which is equivalent to the linearized RG flow. However, in order to determine the 
phase structure and the vortex dynamics of layered systems in a reliable manner one 
has to incorporate the effect of the interlayer coupling which requires corrections to 
the dilute gas result. Our goal here is to show that the the dilute gas approximation, 
i.e. the linearized RG flow, can be used to determine the phase structure of layered 
(multi-flavor) models in the limit of low fugacity (small fermion mass) if the origi- 



2 



nal multi-layer (multi- flavor) model has been rotated in the internal space. For high 
fugacities (large fermion mass) one has to solve the exact RG flow numerically. 



Multi-flavor models. — The bosonized multi-flavor Schwinger model Q can be 
considered as the specific form of a generalized multi-layer sine-Gordon (SG) 
model whose Euclidean action is written as 



S= d\ 



1 1 ^ 

^ ^ n=l 



(3) 



with the 0{N) multiplet = {ipi, . . . ,(Pn). For the specific choice, b"^ = in, and 

(pM^(p^ = iJ? {Yjn=i'^nj , one recovers Eq. Q. The Fourier amplitude related 
to the fermion mass {y ~ m) and the exact relation can be determined by using 
normal-ordering w.r.t. the boson mass. The vacuum angle parameter has to be cho- 
sen as ^ = ±A^7r for i/n > and 6* = for ?/„ < 0. In general, SG-type models 
have also been successfully used to investigate vortex dynamics in 2D or quasi-2D 



superconductors ll24l l25n. Recently, it was shown in Il24t1 that the LSG model with 



a suitable interlayer interaction, 
1 



2- 




^Km-LSG^ =0^[T.^n^n] , (4) 



n=l 



can be used for magnetically coupled layered superconductors where the coupling 
strength between the layers denoted by G and a„ = ±1 are free parameters of the 
model. Based on symmetry considerations [24] any choice with = 1 should 
reproduce exactly the same phase structure, as a consequence, the Fourier ampli- 
tudes (i.e. fugacities) = y for n = 1,2, . . . , N. The frequency 6^ is related 
inversely to the temperature of the corresponding solid-state system. Let us note 
that, different regions of the parameter space have to be considered for the con- 
densed matter and and for the high-energy physics problem. For the description of 
the multi-flavor Schwninger model, one should investigate the phase diagram in the 
two-dimensional plane y — G (for 6^ = 47r) and for the vortex dynamics one has 
to consider the phase structure in terms of the frequency 6^. Let us note, the LSG 
model with magnetic type coupling has a single non-vanishing mass -eigenvalue 
Mjf = NG. Another definition for the mass term of Eq.© 

1 1 



n=l 



is based on the discretization of the the anisotropic 3D-SG model [12611 which has 



been proposed as a suitable model for the vortex dynamics of Josephson coupled 



layered superconductors [27]. Although, it has been shown in [|28l1 . that the LSG 
model with the mass matrix ([5]) cannot be used for Josephson coupled layered su- 
perconductors, in order to distinguish between the two types of mass matrices, in 
this paper we refer to as the Josephson-type interlayer interaction. Let us note 
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that the Josephson-type LSG model can also be considered as a bosonized ver- 
sion of an A^-flavor fermionic model |26, 29, 30fl, but not that of the multi-flavor 
Schwinger model ([T])- In general, the LSG model with Josephson type coupling 
has a single zero and — 1 non-zero mass-eigenvalues, therefore, the Joseph- 
son coupled LSG model is invariant under the particular exchange of the layers 
^ ^N-n+1, hence, ?/„ = Vn-u+i- 



RG approach for multi-flavor models. — In this section we summarize briefly the 
results of the RG analysis of LSG type models discussed in our previous publica- 
tions 1 24, 2^ 28, 29, 30, 31] by means of the approximated form of the Wegner- 
Houghton H differential RG approach (WH-RG). The WH-RG method provides 
us the complete elimination of the modes in the Wilsonian RG method [33] above 
the moving momentum scale k which serves as a sharp cutoff. In principle any types 
of RG methods (see e.g. [134D) can be used to consider the behavior of LSG type 
models. However, the usage of sharp momentum cutoff RG is reasonable since a 
spinodal instability |135i|361] may occur during the flow |l2ji|37i|38|,|39i|40|] and this 
can be used as a signature of spontaneous breakdown of the symmetry ip ^ —ip. 
The WH-RG equation in the local potential approximation (LR^) for the LSG type 
models presented in Refs. \2A, 2^ 2^, 22^ 20^, 21} reads as 



{2 + kdk) 14 M 



1 

47r 



In 



det {dij 



(6) 



where the dimensionless blocked potential Vk = k^'^Vk is introduced and V^^fp) 
denotes the second derivatives of the potential with respect to (pi, ^pj. We make the 
following ansatz for the dimensionless blocked potential of the LSG type models 



N 



(7) 



n=l 



where y„(fc) = k~'^yn{k). Inserting the ansatz ^ into Eq. the right-hand side 
becon ies p eriodic, while the left-hand side contains both periodic and non-periodic 
parts [2a 2^ 30, |3J]. The non-periodic part contains only mass terms, so that 
we obtain a trivial tree-level RG flow equation for the dimensionless mass matrix 

(2 + kdk) M. (k) = 0, which provides the trivial scaling Jk = k~'^J and Gk = 
k^'^G, where the dimensionful interlayer couplings J, G remain constant during the 
blocking. Finally, we recall that in LR\ there is no wave-function renormalization, 
thus the parameter b also remains constant during the blocking. The argument of 
the logarithm in Eq. ([6]) must be positive. If the argument vanishes or if it changes 
sign at a critical value ksi, the WH-RG equation ([6]) loses its validity for k < ksi- 
This is a consequence of the spinodal instability (SI) llssl 361. Below the critical 
scale k < kgi the tree-level blocking relation (see Eq.(13) of 113 ill ) can be used 
to determine the RG flow. In this paper we do not investigate the tree-level RG 
flow of LSG type models but we use the appearance of the spinodal instability as a 
signature of the spontaneous breakdown of the reflection symmetry. 
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In general, the solution of Eq.® can only be obtained numerically, however, var- 
ious approximations of Eq.® are also available in Refs. 11261 |28L |29L |30|, l3l|]. We 
compare two types of approximations, the dilute gas result which is equivalent to 
the linearized RG and the mass-corrected RG flow which incorporates the mass 
term correctly and is able to provide the phase structure of the LSG type models in 
a reliable manner. The linearization of the WH-RG equation ^ in the full potential 
around the UV Gaussian fixed point 114 ill by assuming |5^^Va, | <^ 1, 



N 



[2 + kdk)Vk{v) 



An 



(8) 



n=l 



for the ansatz (|7]) leads to the linearized WH-RG flow Blh exhibiting the solutions 



Vnik) = yn{A) 



A 



(9) 



where y„(A) are the initial (bare) values of the fugacities at the high energy ultra- 
violet (UV) cutoff A. These are the scaling laws valid at the asymptotically large 
UV scales (k ~ A), and being independent of the interlayer coupling (i.e. mass 
terms) predicting a phase structure very similar to that of the massless 2D-SG 



model 
model 




421] . The critical frequency bl = Sir separates the two phases of the 



and the critical temperature is related inversely to the critical frequency 
1/bl BSl . 1371 . kin . Let us note that the linearized WH-RG equations 
obtained in LPA for the = 2-layer LSG model, are the same as those have 



^KTB ~ 



been found in [|27|l in the dilute gas approximation except of the loss of the scale- 
dependence of b due to the usage of the LRA 13 The couplings Jk and Gk are 
always a relevant parameters in the LSG models and, consequently, the lineariza- 
tion, i.e. the asymptotic UV scaling law dH), loses its validity with decreasing scale 
k for any value of b. 



The simplest way to go beyond the linearized (i.e. the dilute gas) approximation 
and to improve the extrapolating power of the UV scaling laws is to take corrections 
into account of the order 0{J/ k"^) for the Josephson and 0{G/ k"^) for the magnetic 
case, which results in the mass-corrected UV scaling laws derived for the LSG 
model in Refs. [24, 2^, 2^. This is achieved by linearizing the WH-RG equation 
in the periodic piece of the blocked potential. 



(2 + fcafc)f/fc(v7i,...,^jv) ~ - 



1 F^{Uk 



47r C 



(10) 



where Uki'-Pi, ■■■■,^n) = Y.n=iyn{k) cos(6v9„) and C and Fi{Uk) stand for the 
constant and linear pieces of the determinant det[5jj + V^*"'] ~ C + Fi(f7fc) +C(f/|). 
Let us first determine the mass-corrected UV scaling laws for the LSG model with 
Josephson type interlayer interaction for = 2. Li this case the solution of Eq. (fTOl) 
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is M,um 




with the initial value y(A) at the UV cutoff k = A. From the extrapolation of 
the UV scaling law Eq. (fTTl) to the IR limit, we can read off the critical values 
bl = IGvr, for = 2. The coupling y is irrelevant for 6^ > bl and relevant for 
6^ < b^. The general expressions for the critical frequency and the corresponding 



critical temperature 12811 read 



bl{N) = SttAT, Tf^sG = ^^= ^KTB^ (12) 
which are determined previously in a similar manner in the framework of the ro- 



tated LSG model in Refs. [|26l. l29Ll30|] . The presence of the coupling J between the 
layers modifies the critical parameter b^ of the Josephson coupled LSG model as 
compared to the massless 2D-SG model. This important modification can only be 
deduced if one goes beyond the linearized (i.e. dilute gas) approximation, e.g. by 
the usage of the mass-corrected UV scaling laws. The similar consideration can be 



done for the LSG model with magnetic type coupling [|24 12911 . Since the layers are 
assumed to be equivalent for the magnetically coupled LSG model, the RG flow 
equations for the fugacities of different layers should be the same {yn{k) = y{k)) 
and the solution can be obtained analytically 



--2 



where y(A) is the initial value for the fugacity at the UV cutoff A and G, b"^ are 
scale-independent parameters. The critical frequency and the corresponding critical 
temperature which separates the two phases of the model can be read directly 

For A^ oo the magnetically coupled LSG behaves like a massless 2D-SG model 
with the critical frequency bl = Stt. In principle, one can try to determine the 
phase structure of the LSG models relying on the dilute gas approximation as it 
has been discussed in Ref. fi43ll for the 2-layer model. However, in this case a 2- 
stages RG procedure is required. In the first step, the real space RG equations are 
integrated out from the UV cutoff (ao ~ 1/A) to the effective screening length 
Ks = 1 / \/2G =1/ V2J, where the topological defects are taken into account with 
full flux. In the second RG step, from Acfr to infinity, an a priori assumption has been 
done by introducing topological excitations with fractional flux and, consequently, 
the predicting power of the RG approach has been weakened. In the next sections 
we show that after an appropriate rotation of layered models in the internal space, 
the dilute gas RG results can be used to determine the phase structure of LSG type 
models without using any a priori assumptions. 
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RG analysis for rotated models. — After performing an 0{N) rotation of the lay- 
ered models which diagonalizes the mass matrix, the rotated models do not have 
interlayer interactions, consequently, the rotated fields can be treated separately. 
Let us note that the rotation has generally been used for coupled two-dimensional 
models, e.g. for the SU(N) Thirring modeljs | and for the 2-flavor 10,0, HQ and 
for the N-flavor Schwinger models IS, 13, 17, 18]. The details of the rotation of the 
A^-layer Josephson coupled LSG model has also been discussed in Refs. ll26l l30n. 
Depending on the number of the non-trivial mass eigenvalues, some of the rotated 
fields have explicit mass terms (massive modes) and the other ones are massless, 
SG-type fields Isi 0, 17]. For example, the dimensionless potential of the rotated 
Josephson coupled LSG model contains — 1 massive fields 



J— rot 



^ 1 ~ 
/ J 2 n n 



+ 



E 

(T1,...,(7]V 



N 

II' 

n=\ 



(15) 



with the rotated 0{N) multiplet = ^ ip^ where Q_ represents the rotation with 
h\ = b'^/N, h"^^^ = IP' / {n{n — 1)) and the integer valued (j„ represent the charges 
of the topological excitations. The rotated magnetic-type LSG model consists of a 
single massive field 



1 



M-rot 



E 

(T1,...,(7JV 



N 
n=l 



(16) 



where = NG and the amplitudes w^^^ ... are different for the Josephson 
and magnetic LSG models. At low energies, below the mass-scale the quantum 
fluctuations are suppressed by the mass terms producing a trivial scaling for the 
massive modes [|4l|], so, the massive modes can be considered perturbatively [SO, 
4 in and they do not influence the phase structure of the rotated models. Therefore, 
one should only consider the remaining massless SG fields in order to determine 
the phases of the layered system. At the lowest order of the perturbation theory, all 
the massive modes are set to be equal to zero. In this case the effective potential for 
the rotated Josephson type LSG model reads as 



Vj_rot(ai) = E' 



i fji h\ ai 



(17) 



and for the magnetically coupled LSG model the effective potential is 



N 



^M-rot(a2, •••,aAr) = E ^f^fT2,...,<7Ar I I 6 



(T2,...,o-]v 



n=2 



(18) 



Let us consider the fundamental modes cr„ = ±1. The linearized WH-RG equation 
^ with V^"" = dlj/k, for the ansatz ([l7]) and leads to the linearized RG flow 
equations 



{2 + kdk)Wi{k) = -^w^{k), {2 + kdk)wM{k) ^ 
4:71 47r 



^ (E WM{k), (19) 
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exhibiting the solutions 

*j(A;)=*j(A)Uj , *M(A;)=*M(A)(^j (20) 

where wj(A) and wmIA) are the initial values for the Fourier amplitudes at the high 
energy UV cutoff A. The critical value of the frequency parameter and the corre- 
sponding critical temperatures are found to be equivalent to Eq.(fT2l) for the Joseph- 
son and Eq. (fT4l) . for the magnetic case. Consequently, the dilute gas RG results for 
the rotated models predict the same layer-dependence of the critical temperature as 
that of obtained by the mass-corrected RG for the original LSG-type models. This 
proves that the dilute gas approximation is suitable to determine the phase structure 
of rotated layered systems for low fugacities. However, for 6^ < 6^, the fugaci- 
ties yn are increasing parameters, consequently, only the exact RG flow is able to 
determine the phase structure of the LSG-type models in a reliable manner. 

Exact RG flow for the N = 2 flavor model. — Since the results (fT2l and [141) have 
been established by an UV approximated RG method, it is certainly worthwhile to 
confirm the analysis by a numerical calculation of the exact RG flow. Moreover, if 
one considers the appearance or non-appearance of spinodal instability during the 
blocking which can be used as a signature of spontaneous symmetry breakdown of 
the reflection symmetry, the numerical solution of the exact WH-RG equation is 
required. Indeed, the full RG analysis of the 1 -layer LSG model at 6^ = An dis- 



cussed in Ref. i21|] provides us the tool to investigate the symmetric and the sym- 



metry broken phases of the 1 -flavor massive Schwinger model. The critical value 
of the ratio i^)c = 0.311 which separates the two phases of the 1 -flavor model has 
been determined by the exact RG method which coincides with the results of other 
calculations (see e.g. [10]). If g ^ m, i.e. below the critical ratio, the spinodal 
instability does not appear during the RG flow, therefore, the reflection symmetry 
remains unbroken. This is the consequence of the trivial scaling of the Fourier am- 
plitudes below the mass-scale. However, for g <^ m, i.e. above the critical value, 
the spinodal instability always appears. One may assume a similar phase structure 
for the multi-layer model, however, it has been argued in the literature [si [isi \^\ 
that the low-energy behavior of the multi-flavor Schwinger model is different de- 
pending on whether = 1 or > 2. Our aim here is to clarify this issue by 
the numerical solution of the exact WH-RG equation derived for the 2-layer LSG 
model. We will show that spinodal instability always appears for the 2-layer LSG 
model for b"^ = A-n. 

We determine numerically the dimensionless effective potential Kff(<^i) '^2) for the 
double-layer LSG model as the limit A; ^ of the dimensionless blocked potential 

Vki'^i, f2) = ^Gk{ip2 - <^lf + Ukifl, (P2) , (21) 
where Uki^i, ip2) is an arbitrary periodic function of the fields (with Z2 symmetry) 
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including all the Fourier modes generated during the RG flow. In order to consider 
the effect of the higher Fourier modes, which were not taken into account in the 
previously utilized linearized and mass-corrected linearized WH-RG approach, we 
use the following ansatz for periodic part of the blocked potential 



Uk = uoi{k)[cos{b(pi) + cos{b(p2)] + uii{k) cos{bipi) cos{bip2 
+ viiik) sin(6v9i) sin(6v92), 



(22) 



with the fundamental mode mqi = mio = V- Inserting the ansatz (I2TI) into the exact 
WH-RG equation Q and separating the periodic and non-periodic parts, one ar- 
rives at the RG equation for the periodic part, see Eq.(17) of Ref.f3lll. For technical 
reasons, it is more convenient to consider derivative of the WH-RG equation with 
respect to one of the field variables. By Fourier decomposition, this RG equation 
can be reduced to a set of ordinary differential equations for the couplings uqi, uu 
and Vii, 



A 



where Di 



DkUii 

\DkVii ) 
kdk 



( 



47r 



-2(1 + G)uqx + b^UoiUn 



\ 



-2(1 



G)un + 2Gvn + b^u^ 



11 



V 



(23) 



-2{1 + G)vn + 2Guu J 



and element of the matrix A are A 



4i3 = ^31 = T^oi^ii, and = ^32 



-2(1 + 2G) + 

2G), 42 = A21 = + G)uoi - yMoiMii, 
0. We invert the matrix A and solve the 
RG flow equations for the Fourier amplitudes uqi, uu and vu numerically, by a 
fourth order Runge-Kutta method, whose numerical stability was verified by vary- 
ing the step size. The main advantage of the numerical solution of the WH-RG 
flow equation (l23l) is that all the non-linear terms are kept. The numerically de- 
termined scaling of the couplings uoi{k), un{k) and vu{k) can be compared to 
the corresponding approximate UV scaling laws. The results are the foUowings. 
The mass-corrected UV scaling law (fTOl) for the fundamental mode (mqi) coincides 
with the numerically obtained one. For the example, for 6^ = 127r, the deviation 
between the numerical solution of the exact WH-RG equation and the solution of 
the mass-corrected UV linearized WH-RG equation at the scale A; = 1.0 x 10^^ 
is 9.69 X 10^^. This coincidence demonstrates that the flow of the fundamental 
coupling Uqi is well described by the mass-corrected UV scaling law if no spinodal 
instability appears during the blocking. Therefore, we find that for low fugacities 
(small fermion mass), the phase structure of the double-layer LSG model obtained 
numerically is the same as that predicted by the extrapolation of the mass-corrected 
UV scaling laws, because the higher harmonics do not modify the scaling of the 
fundamental mode (see Figs. [T] and [21). The IR behavior of the higher harmonics 
and consequently the effective potential can only be determined by the numerical 
solution of the exact WH-RG equation when all the non-linear terms are kept. For 
example, the UV approximated RG flow is not able to determine the sign changes 
(see the peaks of the dotted and dashes lines in the figures) of the higher harmonics. 
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Fig. 1. The exact RG scaling of the dimensionless coupling constants uqi, uu and of the 
double-layer LSG model is represented graphically for various frequency parameters. The 
mass scale is M/A = VO.002. The full, dotted and dashed lines correspond to uqi, uh and 
vu respectively. Spinodal instability (/csi/A) appears below the mass scale. For increasing 
value of the momentum scale A;si tends to zero and vanishes at 6^ = 167r. Belowfei the 
WH-RG equation (O looses its validity and the tree-level RG relation (Eq.(13) of 113 IH ) has 
to be used which is not discussed here. 
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Fig. 2. The scaling of the dimensionless coupling constants of the double-layer LSG model 
is shown for two different frequencies = IOtt and b^ = ISvr and for the mass scale 
M/A = VO.002. The UV (A/ < k) and IR (A; < M) regions are separated by the mass 
eigenvalue. The scaling of the fundamental mode ttQi (full line) depends on the frequency, 
increasing (relevant) for 6^ < IGn and decreasing (irrelevant) for b^ > IGvr. The initial 
values and consequently the UV scalings are different for uu (dotted) and vu (dashed) but 
in the IR region the trajectories are coincides resulting in a trivial tree-level scaling (~ k^^) 
which is independent of the frequency. 

According to the exact RG flow, in the IR limit uu and vu coincide independently 
of their initial values at the UV cutoff (see Fig. O. Therefore, the low-energy ef- 
fective potential of the 2-layer LSG model determined by the exact RG approach 
is 



off 



+ uoi[cos{b(pi) + cos(6v92)] + wii cos{bipi - b(p2), (24) 



with wu = ^tii = Vu which corresponds to the massive mode in the rotated form 
of the model, hence, Wu has a trivial tree-level IR scaling (~ which is inde- 
pendent of the frequency 6^. The numerical solution of the exact WH-RG equation 
shows (Fig.[T]) that for the 2-layer LSG model with 6^ = Att, the spinodal instability 
(SI) always appears during the RG flow in contrary to the 1 -layer model where the 
appearance of SI can be avoided for sufficiently small initial value for the fugacity. 
For the 1 -flavor model for low fugacity, the SI cannot be detected below the mass 



10 



scale since the presence of the mass term predicts a trivial tree-level scaling for 
all the Fourier amplitudes. For the 2-flavor model, below the mass scale, only the 
massive modes have trivial scalings but not for the fundamental one (mqi). Hence, 
the appearance of SI is unavoidable. 



Summary. — In this paper we investigated the phase structure and the low-energy 
behavior of the bosonized multi-flavor Schwinger model ^ by means of the 
Wegner-Houghton RG method. The Bose form of the multi-flavor model consists 
of 2D sine-Gordon fields coupled by an appropriate mass matrix (HI) which has 
also been used to describe the vortex properties of magnetically coupled layered 
superconductors ll24ll . Another definition (see Eq.®) for the mass matrix of the 
multi-flavor i.e. layered sine-Gordon model has also been discussed which is based 
on discretization of the 3D sine-Gordon model [26]. Let us note, that the layered 
model with the mass matrix ^ can in principle also be considered as a Bose form 



of a fermionic model as it has been argued in [|26l . 13011 . however not that of the 
multi-flavor Schwinger model 



It has been shown that in the limit of small fermion mass the linearized RG flow (i.e. 
the dilute gas approximation) is sufficient to determine the phase structure of the 
multi-flavor model in a reliable manner, if it has been rotated by a suitable rotation 
in the internal space which diagonalizes the mass matrix. This receives important 
application in condensed matter physics where the usual RG techniques are based 
on the dilute gas approximation. For example, using the idea of rotation the vortex 
dynamics of magnetically coupled layered superconductors can be considered by 
means of dilute gas RG methods and no two-stages RG ll43n is required. The cal- 
culation of the exact Wegner-Houghton RG flow of the 2-flavor (i.e. 2-layer) lay- 
ered sine-Gordon model, by a numerical approach including higher-order Fourier 
modes confirms that the mass-corrected UV scaling law (i.e. the linearized RG flow 
for the rotated model) is sufficient to determine the phase structure of the layered 
sine-Gordon model in the low fugacity (small fermion mass) limit. 



The rigorous RG study of the layered sine-Gordon model verifies that the low 
energy effective potential of the bosonized multi-flavor Schwinger model is a sine- 
Gordon type model which undergoes a KTB-type phase transition at the flavor- 
number dependent critical frequency h1{N) = SttN/ {N — 1). In the limit — oo, 
the layered sine-Gordon model tends to the two-dimensional sine-Gordon theory 
with the critical frequency = 8n. Therefore, in the large limit, the low energy 
behavior of the bosonized multi-flavor Schwinger model becomes independent of 
the boson mass i.e. the coupling g. This is consistent with the flavor-dependence 
of the chiral condensate < i/j^ijj >~ 77i(^"i)/(^+i)(j(2/(^+i)^ and the mass gap 
Mgap ~ mN/(N+i)gi/{N+i) ^Yiich are independent of ^ if A^ ^ oo [17]. The nu- 
merical solution of the exact RG equation shows that in case of the 2-flavor (i.e. 
2-layer) layered sine-Gordon model with 6^ = An, the spinodal instability always 
appear during the RG flow in contrary to the 1 -layer model where the appearance 
of spinodal instability can be avoided for sufficiently small initial value for the fu- 
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gacity. Consequently, for > 1 the reflection symmetry always suffers breakdown 
in both the weak and strong coupling regimes, in contrary to the = 1 case, where 
it remains unbroken in the strong coupling phase. 

Finally, let us mention that the extension of the RG analysis presented in this paper 
can also be used to consider the 6'-dependence of the multi-flavor model and to map 
out the phase structure of the multi-frequency sine-Gordon type models [i44ll . 
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